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Abstract 



The recent demonstration of current singularity formation by Low et al. ass umes 
that potential fields will remai n potential under simple expansion or compression (|Low 
20061 . 120071 ; iJanse &; Lowll2009l ). An explicit counterexample to their key assumption is 
constructed. Our findings suggest that their results may need to be reconsidered. 



1. Introduction 

The theory of force-free magnetic fields is of great interest in many astrophysical situations 
(jMarshl Il996l : ISturrockl Il994l ). For instance, in the quiescent solar corona, the magnetic energy 
dominates the thermal energy, the kinetic energy, and the gravitational potential energy. Under 
these conditions, the Lorentz force J x B of the magnetic field is approximately zero, therefore the 
field is force- free to the lowest order. A force-free field satisfies 



J = V x B = AB; 



(1) 



that is, the electric current J is parallel to the magnetic field B everywhere. From the conditionV • 
B = 0, we obtain B • VA = 0. 

The magnetic field lines in the solar coronal loops are deeply anchored in the dense photosphere, 
which is continually moving due to the convective motions in the solar interior. As a result, the 
coronal magnetic field is constantly changing in response to the photospheric motions. Since the 
typical turnover time of photospheric motions (^hundreds of seconds) is much longer than the 
Alfven transit time (~tens of seconds), the solar coronal loops remain approximately force-free for 
all time, as long as they are quiescent. In other words, to a good approximation, the photospheric 
motions carry the coronal magnetic field through a sequence of quasi-static, force-free equilibria. 
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It is, of course, possible that the force-free field may become unstable to magnetohydrodynamic 
(MHD) instabilities. When instabilities occur, the coronal loops will no longer remain quasi-static, 
and subsequently some violent events, such as coronal mass ejections (CMEs), may ensue. 

The solar corona is highly conducting. To a very good approximation, the plasma resistivity is 
negligible. One of the most important problems in solar physics is to explain how the solar coron a 



could be heated to millions of degrees when the resistive dissipation seems so weak ()Klimchukll2006l ) . 



( Parker 


1972, 


1983. 


1979, 


1994) 



ideal plasma has a natural tendency to form tangential discontinuities, that is, current singularities, 
when the field line topology is rendered sufficiently complicated by the boundary motions of the 
photosphere. In reality, because of the finite resistivity, the ideal current singularities will be 
smoothed out and become thin current filaments. The current filaments are so intense that even 
the tiny resistivity of solar corona can lead to significant dissipation. Parker suggests that the 
dissipation from ubiquitous thin current filaments could be the energy source of coronal heating. 



Janse fc Low 



20091 ). Low and coauthors try to demonstrate unambiguously the formation of current 



Parker's claim has stimulated considerable debate that continues to this day (van Balleeooiien 


1985 




Zweibel & Li 


1987; Antiochos 1987; LonecoDe & Strauss 1994; Ne & Bhattachariee! 


1998; 


Craie: & Sneyd 


2005). without apparent consensus. Recentlv. in a series of papers ( 


Low 


2006, 


2007; 



singularities using what they call "topologically untwisted fields", which, in the context of force- 
free field, is synonymous with potential fields. That is, they limit themselves to a special subset 
of force-free fields with A = in Eq. (pQ), which correspond to current-free fields. In this case, 
the magnetic field can be expressed as B = V% with some potential Xi where x satisfies V 2 x = 
because B is diver genceless. In a simply connected, compact domain, x 1S uniquely determined 
up to an additive constant if the normal derivative of x 1S given on the boundary. Therefore, a 
potential field is uniquely determined by prescr ibing the normal com ponent of B (= n ■ Vx, where 
n is the unit normal vector) on the boundary. iJanse fc Low! (J2009J) consider a potential field in a 
cylinder of finite length as an initial condition; the normal component of the field is nonzero only 
on the top and the bottom of the cylinder. The cylinder is then compressed to a shorter length. 
Because the whole process is governed by the ideal induction equation, the normal component 
of the field remains the same. By making a key assumption that the field will remain potential 
during the process, they circumvent the complicated problem of solving for the various stages of 
quasi-static evolution, and simply calculate the magnetic field in the final state as the potential 
field which satisfies the boundary conditions. They then numerically calculate the magnetic field 
line mapping from one end to the other, for both the initial and the final states. They find that for 
a three-dimensional (3D) field, the field line mapping, and hence, the topology, is changed. They 
claim that current singularities must form during the process; furthermore, singularities may form 
densely due to the ubiquitous change of the field line mapping. 



The assumption that the field will remain potential during the compression, however, is not 
proven in their work. Intuitively it looks plausible, as parallel current is related to the relative twist 
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between neighboring field lines. It would appear to require a nonzero vorticity on the boundary, 
which is absent in a simple expansion or compression, to twist up the field. However, parallel current 
is related to the relative twist between neighboring field lines only in an average sense; the relative 
twist between a specific pair of field lines is actually independent of the parallel current. Although 
evolutions following the ideal induction equation preserve the footpoint mapping of each individual 
field line, it does not follow that parallel current will not be generated during a simple expansion 
or compression. In this paper an explicit counterexample is constructed in a two-dimensional (2D) 
slab geometry; therefore the demonstration of current singularity formation by Low et al. is put in 
doubt. 



2. Two-dimensional Counterexample 

Consider a 2D configuration in Cartesian geometry, with x the direction of symmetry. A 
general magnetic field may be expressed as 

B = B x x + x x W, (2) 

where B x = B x (y,z) and ip = ip(y,z). For a force-free field, B x = B x (ip) and ip must satisfy the 
Grad-Shafranov equation: 

9 „ dB r , . 

W— fl.-^. (3) 

The A for a force-free field in Eq. ([1]) is equal to —dB x /dip in this representation. If the field is 
potential, then B x = const and V 2 ip = 0. 

Let us consider a force- free field enclosed by two conducting boundaries at z = and z = L, 
with all field lines connecting one end plate to the other. The footpoint mapping from one end to 
the other is characterized by ip on the boundaries, as well as the axial displacement along x when 
following a field line from the bottom to the top: 

Here ds is the line element on the y — z plane, and the subscript ip indicates that the integration 
is done along a constant ip contour. Notice that h{ip) = d$>/dip, with $>(ip) the axial magnetic flux 
function up to an additive constant. The displacement function h(ip) resembles the safety factor 
q(ip) in toroidal geometry. 

Let the initial condition be {B X Q,ipo) with a system length Lq. Suppose the system has un- 
dergone a simple expansion or compression such that the system length changes to some L ^ Lq. 
Since ideal evolution preserves the footpoint mapping, the flux function i/j on the boundaries and 
the axial displacement h(ip) must remain unchanged. To determine the final force-free equilibrium, 
we have to solve the two coupled equations 

(5) 
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subject to 7p(y,0) = tpo(y,0), ip(y,L) = tpo(y,Lo), where the axial displacement h(ijj) is determine 
from the initial condition. The set of coupled equations (El) and (El) is called a generalized differential 



equation, which in general requires numerical solutions ( Grad et al.lll975l ). 



As a simple example, consider the following initial condition: ipo = y + esin(y) cosh(z — 1/2), 
BxO — 1) an d Lq = 1, where e <C 1 is a small parameter. This initial field is potential. Now 
the question is, will the final field remain potential when the system is expanded or compressed 
to another length L? We may first assume that it will, and see if that leads to a contradiction. 
Assuming that, we have to solve V 2 ip = subject to boundary conditions ip(y,0) = ip(y,L) = 
y + esin(y) cosh(l/2), and B x = const. The solution for ip is 

tp = y + easin(y) cosh(z — L/2), (7) 

where 

_ cosh(l/2) 

" cosh(L/2) ' l8j 

The constant B x may be determined by the conservation of axial flux, which yields B x = 1/L. To 
be consistent with the initial footpoint mapping, the final state (ip,B x ) must give the same axial 
displacement as (ipo, B x q). We examine this along the field line y = 0, corresponding to tpo = tp = 0. 
For the initial field, 

1 1 

1/2 



h(ip = 0) = / — — — — dz 

1 + ecosh(z — 1/2) 



and for the final field, 



1 



-1/2 1 + e cosh z 



L 



7 dz>; (9) 



h(ip = 0) = - / — — -dz 

L J 1 + ea cosn(z — L/2) 



1/2 1 

1 TTt — 7\d z> ■ (10) 

-1/2 1 + eacosa(Lz ) 

The condition that the axial displacement remains unchanged implies 

' 1 / 2 acosh(Lz') — coshz' , 

7 TTT^ 7 T^dz = 0- (11) 

1/2 (1 + ecoshz') (1 + eacosh(Lz ')) 

However, it may be shown that the integrand of Eq. pip for z' £ (—1/2,1/2) is positive when 
L < 1 and negative when L > 1; therefore the condition pip cannot be satisfied except for the 
trivial case L = Lq = 1 and we have a contradiction. 
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Having shown that a potential field cannot be the final force-free equilibrium, the next question 
is whether a smooth, non-potential equilibrium exists. To demonstrate the existence of a smooth 
solution, let us consider a perturbative solution 

V> = y + eVi+0(e 2 ), (12) 

with the following ansatz 

ip! = sin(y)g(z). (13) 

From the initial condition, we have 

h(ip) = 1 - 2esinh^cos(^) + 0(e 2 ). (14) 

And Eq. ©, to 0(e), yields 



BM = hty) ( I" L ' T^-rrds 



=o |V^| 

h(if)) ^L-ecos^J g(z)dz + 0(e 2 ) 



1 e cos ip / 1 



L L 

The right hand side (RHS) of Eq. (J5J is 

dB x 



i ^ L s(*)<fe - 2 sinh + 0(e 2 ). (15) 



± £ g(z)dz- 2 sinh^J + 0(e 2 ). (16) 



Hence, Eq. ©, to 0(e), is 

</'(z)-g(z) = ^ (jj\{z)dz - 2 sinh . (17) 

And the boundary condition for ^ requires g(0) = g(L) = cosh(l/2). The solution can be readily 
found as 

V 



g(z) = ci cosh Iz - —J + c 2 , (18) 
where 

L (L 2 + l) cosh(l/2) — 2Lsinh(l/2) 

Cl = Z (L 2 + 1) cosh(L/2) - 2 sinh(L/2) ' (19) 
and i L 

C2 = cosh — — ci cosh — . (20) 
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Therefore, we have shown that a smooth solution could indeed be found, at least to 0(e). From 
Eq. (|15p . the axial magnetic field is 

BM = (|ci sinh | + c 2 - 2 sinh ±) + 0(e 2 ). (21) 

And finally, 

dB x esin-0 / 2 , L , 1\ ^. 9 . . . 

A = "d# = (l Cl Smh 2 + C2 " 2 2 J + 0(6 )■ (22) 

Since A 7^ almost everywhere when L 7^ 1, a smooth electric current is generated through out the 
whole domain. 

We have checked the perturbative solution by the numerical solution of the generalized dif- 
ferential equation. We find good agreement for small e. For larger e the numerical solution does 
not quite agree with the solution, as expected. Furthermore, we have been able to find smooth 
numerical solutions for a wide range of L and e. Figure [T] shows LB x (ip) as a function of ip for 
two representative lengths L = 0.1 and L = 0.9, with e = 0.1 and e = 0.5. Solid lines are numer- 
ical solutions and dashed line are perturbative solutions. The e = 0.1 case shows good agreement 
between the perturbative and numerical solutions. For e = 0.5 the agreement is only qualitative. 

3. Summary and Discussion 

In summary, we have constructed an explicit counterexample to the assumption that potential 
fields remain potential under simple expansion or compression. For our 2D system, the reason 
that in general the field will not remain potential is quite simple. A potential field has to satisfy 
V 2 ijj = 0, which uniquely determines for the given boundary value of ip. This ijj, however, has to 
satisfy an extra constraint to be self-consistent. That is, when substituting ip into the RHS of ([6|), 
we should get a constant. This is a very stringent constraint and in general will not be satisfied, 
as we have demonstrated with a simple example. 

It should be pointed out that the domain under consideration is not compact, therefore spec- 
ifying the normal component of the magnetic field on the boundary does not uniquely determine 
the potential field. One can, for example, add a constant B y or B x to the field without changing 
the normal component of B. In this regard, our configuration does not belong to the same class 
of configurations as considered by Low et al. However, the purpose of this work is to demonstrate 
that a potential field can evolve into a non-potential field by simply changing the system length, 
therefore it should not be taken for granted that the final field has to remain potential. Since this 
is the key assumption in the demonstration of current singularity formation by Low et al., our 
findings suggest that their conclusion needs to be revisited. Furthermore, Eq. (|22p shows that a 
parallel current is generated almost everywhere. This may account for why Low and coauthors find 
ubiquitous change in the footpoint mapping. Although our analytic demonstration of the existence 
of a smooth solution is based on perturbation theory, the demonstration is supported by numerical 
solution even when perturbation theory is not strictly valid. 
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The question posed by title "do potential fields develop current sheets under simple compression 
or expansion?" remains open at this point. Indeed, we cannot preclude the possibility that current 
sheets may form during the process. However, even if that is the case, a smooth current density 
is likely to precede the formation of current sheets. This raises the overall question as whether 
limiting to potential fields is a viable option in solving Parker's problem. The answer is probably 
no. 
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Fig. 1. — LB x (ip) as a function of ip for two representative length L = 0.1 and L = 0.9. Panel (a) 
shows the case e = 0.1 and panel (b) shows the case e = 0.5. Solid lines are numerical solutions 
and dashed line are perturbative solutions. The e = 0.1 case shows good agreement between the 
perturbative and numerical solutions. For e = 0.5 the agreement is only qualitative. Notice that 
B x = 1 when L = 1. Since B x {ip) becomes non-constant for L / 1, the field becomes non-potential. 



